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We investigate the azimuthal asymmetry v2 of produced pions in pp and pA collisions at both 
RHIC and SPS energies. In our approach, based on the pQCD parton model and the light-cone 
QCD-dipole formalism, the azimuthal asymmetry results from a correlation between the color-dipole 
orientation and the impact parameter of the collision. We introduce the color-dipole orientation 
within an improved Born approximation and the saturation model which satisfies available DIS 
data, showing that the azimuthal asymmetry of partons and pions is very sensitive to the choice 
of the model, and that it is reduced in the saturation model. We find that v2 of quarks and 
gluons in parton-nucleus collisions have very different patterns. The azimuthal asymmetry of gluons 
. in gluon-nucleus collisions can be negative at small transverse momentum, changes the sign and 

' becomes positive at high transverse momentum. The azimuthal asymmetry of quarks in quark- 

I nucleus collisions is positive at all values of transverse momentum. We find that the azimuthal 

anisotropy v2 of produced pions in both pp and pA collisions is positive, albeit rather small. 
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I. INTRODUCTION 



'" ^ One of the most important discoveries in uhrarclativistic nucleus-nucleus collisions at the Relativistic Heavy Ion 

I Collider (RHIC) was the observation of a large elliptic flow which is as large as the one predicted by ideal-fluid 
Qh, models (which assume zero viscosity) 0, 0] . Microscopically, the elliptic flow results from the interactions between 
^ ■ the produced particles, and therefore carries information about the dense matter produced at RHIC. 

Various theoretical calculations [2] support the notion that collective flow is perhaps generated early in the nucleus- 
nucleus collision, and is present at the partonic level, before partons coalesce or the hadronic fragmentation stage. 
On the other hand, it has been shown that two-body interactions between partons cannot by themselves generate 
sufficient flow to explain the observations, unless partonic cross sections are artificially enhanced by more than an 
' order of magnitude over perturbative QCD predictions (see, however, the recent development on this line in 
\ Ref. @). This indicates that the quark-gluon matter created at RHIC is strongly interacting, unlike the type of 
. weakly interacting quark-gluon plasma expected to occur at very high temperatures on the basis of asymptotic 
freedom. Ultrarelativistic heavy ion collisions probe QCD in different regimes at different stages of the collisions. To 
be able to understand the underlying dynamics of the collisions and in particular the observed elliptic flow, one has 
' to learn how to disentangle between the contributions of different physical subprocesses and their contributions to the 
[ observed azimuthal asymmetry. 
|!^ ' We have proposed recently a scenario which produces an azimuthal asymmetry coming from the initial stage of 
relativistic nuclear collisions '^^^ approach, the azimuthal asymmetry is related to the sensitivity of parton 

multiple interactions to the steep variation of the nuclear density at the edge of the nuclei via the color dipole 
^ i orientation. This effect is also present in elementary pp reactions where the correlation between the color-dipole 
orientation and the impact parameter of the collision leads to an azimuthal asymmetry. We have recently computed 
the azimuthal asymmetry of the prompt photons in nuclear collisions coming from this mechanism Q. 

In this paper we use this idea to calculate the azimuthal asymmetry of pions produced in pp and pA collisions. 
Measuring V2 of the produced particles in pp and pA collisions is a challenge for experimentalists, partly due to 
the difficulties associated with the identification of the reaction plane. Nevertheless, measurements of the azimuthal 
correlations from pp, dAu and AuAu collisions at RHIC indicate that the azimuthal asymmetries in pp and pA 
collisions can be quite different from those in AA collisions [H . 

Here we systematically study the azimuthal asymmetry of quarks and radiated gluons in parton-nucleus collisions, 
and show how they contribute to the azimuthal asymmetry of produced hadrons in pp and pA collisions. The basic 
tools of our calculations include the pQCD parton model combined with the light-cone QCD-dipole formalism. Such 
a study helps us to understand how the azimuthal asymmetry of the produced particles evolves from elementary pp 
collisions to cold nuclear matter in pA collisions and then to a hot quark-gluon plasma in AA collisions. The azimuthal 
asymmetry of the produced particles in peripheral AA collisions or at high pT is expected to be similar to those in 
pA and pp reactions. Therefore, this study may be also relevant for the physics of AA reactions and can be used as a 
baseline for jet-quenching models, since in pp and pA collisions no hot and dense medium is created. 

This paper is organized as follows: in Sec. I and II we introduce the color-dipole orientation in an improved Born 
approximation 11] and within the saturation model of Golec-Biernat and Wiisthoff [12]. In Sec. IV, we calculate the 



2 



azimuthal asymmetry of quarks in quark-nucleus collisions. In this section, the broadening of a projectile parton in 
pA collisions will be introduced. This is the key ingredient of the Cronin effect and hadron production in pA collisions. 
In Sec. V, we study the azimuthal asymmetry of the radiated gluons in gluon-nucleon and gluon-nucleus collisions. 
In Sec. VI, we discuss hadron production in the short- and long-coherence regimes. We introduce two very different 
schemes for hadron production in pp and pA collisions: (1) pQCD-improved factorization and (2) the light-cone QCD- 
dipole factorization scheme. The numerical results and discussions are given in Sec. VII. Some concluding remarks 
are given in Sec. VIII. In an Appendix, we give some details of the calculations carried out in Sec V. 



II. COLOR DIPOLE ORIENTATION IN BORN APPROXIMATION 



A colorless qq dipole is able to interact only due to the difference between the impact parameters of the q and q 
relative to the scattering center. Therefore, a qq fluctuation cannot be produced if both q and q have the same impact 
parameter b from the target, even if their transverse separation r is not zero. In terms of the partial elastic amplitude 
fqq{b, r), it means that the vectors r and b are correlated. This can be seen in a simple example of a dipole interacting 
with a quark in Born approximation. The partial elastic amplitude up to a factor M reads. 
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where Kf){x) is the modified Bessel function and we introduced an effective gluon mass nig to take into account some 
nonperturbative effects. It is obvious from the above expression that the partial elastic dipole amplitude exposes a 
correlation between r and 6, and the amplitude vanishes when 6 • r = 0. 

The Born amplitude is unrealistic, since it leads to an energy independent dipole cross section aqq{r,x) where 
X corresponds to the Bjorken variable in DIS. The energy dependence of the dipole cross section is generated by 
additional radiation of gluons, which can be resummed in the leading ln{l/x) approximation. As a first step to 
improve the above approximation, we first obtain the dipole cross section from Eq. ([1]) at small r [llj . 



lim (To 
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the unknown coefficient J\f can be then obtained by comparing the above expression with the small r expansion of the 
Golec-Biernat and Wiisthoff (GBW) dipole cross section Notice that at small r there is no logarithmic term in 
the GBW dipole cross section, therefore we fix the mean value of f, which depends on the process under consideration. 
Thus, we obtain 



27r (0.62-ln(mgf))i?2(a;)^ 



(3) 



where ctq = 23.03 mb, Ro{x) = 0.4 fm x {x/xo)°-^'^'^ with xq = 3.04 x 10-"*. Notice that the amplitude given in Eq. ((H) 
is for the qq dipole colliding with a quark target, not with a nucleon. Although by fixing the coefficient J\f with the 
GBW cross section some of the missing effects are incorporated, still the above approximation is rather crude. In the 
next section, we introduce the dipole orientation without using any approximation. 



III. COLOR DIPOLE ORIENTATION IN THE SATURATION MODEL 



Here, we introduce the color dipole orientation within the phenomenological saturation model of GBW, which 
includes contributions from higher order perturbative corrections as well as non-perturbative effects contained in DIS 
data. The dipole elastic amplitude of a qq dipole colliding with a proton at impact parameter b is given by Q 

where we defined = \/cts{q'^)as{Q''^), and T{x, q, q') is the generalized unintegrated gluon density. The fractional 
light-cone momenta of the quark and antiquark are denoted by /3 and 1 — respectively. In the 2g model, Eq. ([1]), 
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we assumed that that q and q have equal longitudinal momenta, i.e. they are equally distant from the dipole center 
of gravity, which corresponds to the case with a parameter (3 = 1/2. The generalized unintegrated gluon density was 
proposed in Ref. Q assuming that the transverse momentum distributions of the two gluons do not correlate, except 
for the Pomeron-proton vertex, which is function of the total momentum transfer, and it is given by 
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This generalized unintegrated gluon density is related to the diagonal one by J-{x, q,q' — q) = J-'{x, q). We assume here 
that the transverse momentum dependence of the dipole-proton elastic amplitude has a Gaussian form. Comparison 
with the saturated form of the dipole-proton cross section <Tqq{r, x), 



= 2 / d%lmf^Jb,f,x,P) 



(6) 
(7) 



fixes the form of J^{x, q, q') and the values of the parameters, except for Rj^(x). 

Knowing the generalized unintegrated gluon density, we can now perform the integration in Eq. (j4]) and obtain the 
partial elastic dipole-proton amplitude. 
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with the notation B{x) — R%{x) + Rl{x)/8. 

To fix the function R%(x) we use another observable, the t-slope B'^i'{x,r) of the elastic dipole-proton cross section 
(at f = 0) in the limit of vanishingly small dipole, r 0. In this limit. 



Bl1{x,r^Q) = B{x) + -Rl{x). 



d^bb'^lmf^Jb,f,x, (3 = 1/2) =B{x) + 
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In this expression we fixed (3=1/2 for the sake of simplicity. Equation ^ can be compared with the slope of the cross 
section of elastic electroproduction of p-mesons measured at HERA at small x and high . It was observed that at 
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10[ , which can be compared with 
Therefore, 



Q'^ ^ 1 GeV^ the slope saturates at the value Bjtp^pp{x, ^ 1 GeV"^ 

our result Eq. ^ in the limit r ^ 0, since at high the effective size of the dipole is vanishingly small, 
we have R%{x) = B^*p^pp{x, > 1 GeV^) - \RI{x). 

Notice that the expression Eq. ([5]) at r ^ exposes the property of color transparency [ll|j f^q(b,r,x, (3) oc r^. 
It also goes beyond the usual assumption that the dipole cross section is independent of the light-cone momentum 
sharing (i. Although, the partial amplitude Eq. ^ does depend on (3, this dependence disappears after integration 
over impact parameter b as shown in Eq. ([7]). From Eq. ([8]), it is seen that when the transverse dipole size r and the 
impact parameter b become comparable in size then the orientation becomes important. For very small r or b. the 
dipole orientation is not present. The partial dipole amplitude behaviour also changes with the parameter j3 

One should note that the GBW model is a simple parametrization which has some restrictions. In particular, the 
model exhibits no power-law tails in momentum space in contradiction with QCD. Besides, it does not match the 
QCD evolution (DGLAP) at large values of . Therefore, one should be cautious applying this model at very high 
transverse momenta accessible at the energies of LHC. 

In the above we relied on the saturation GBW model, which depends on Bjorken x. However, for soft reactions the 
cm. energy s, rather than Bjorken x, is the proper variable. Similar to the GBW model, the s-dependent dipole cross 
section with a saturated shape fitted to data on DIS at not high, and to real photo-absorption and photoproduction 
of vector mesons, was introduced in Ref. [3, 



j(r, s) CTo(s) 
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where Ro{s) = 0.88 fm (so/s)"'"'^'' with Sq = 1000 GeV^. The normahzation factor ao{s) is fixed by demanding that the 
pion-proton total cross section be reproduced, that is / (Pr |\['^(r)ptT^q(r, s) = (T^JJ(s), where the pion wave fmiction 
squared integrated over longitudinal quark momenta has the form, 

with a mean pion charge radius squared {rl^f^-K = 0.44 fm^ 15]. In this way, the normalization factor (Tq{s) is 
determined, 

We employ the parametrization of the fit in Ref. [16] for the Pomeron part of the cross section (t^q'J(s) = 
23.6(s/so)°"^ mb, where sq = 1000 GeV^ 

We assume that for soft processes one can switch from x- to s-dependence, keeping the same functional form of 
the dipole amplitude Eq. ([5]) but adjusting the parameters and cro('S) to observables in soft reactions. The 

first condition is that the s-dependent dipole partial amplitude reproduces the s-dependent pion-proton cross section. 
Another condition is the reproduction of the slope at < = 0, B^f{s) = \{b'^). These conditions will be satisfied by the 
following replacements: 

lmfg{b,r,x,P) ^ Im/^^(6,f,s,/3), 

Ra{x) =^ i?o(s) = 0.88fm(so/s)°-" 

ao ao{s) = aZis) (l + fpj;) ' (KST model) (13) 

where we use a Regge parametrization for the elastic slope, B^f{s) = Bq + 2aip ln(s//Lt^), with Bq = & GeV^^, 
ajp = 0.25 GeV~^, and ^J? — \ GeV^. In what follows, we call the s-dependent dipole amplitude KST model. 



IV. V2 OF QUARKS 

Multiple interactions of projectile quarks in the target may proceed coherently or incoherently. In the former case 
the multiple interaction amplitude is a convolution of single scattering amplitudes, and in the latter case one should 
convolute differential cross sections, rather than amplitudes. The condition of coherence is exactly the same as in 
classical optics, namely the maximal longitudinal distance between different scattering centers should not considerably 
exceed the so called coherence length. 



Ic = (14) 
Pt 



where Eq is the quark energy in the nuclear rest frame and pr is the total transverse momentum of the quark 
accumulated from the multiple rescattcrings. Notice that at mid rapidities we have Eq = pT^fs/mM- Thus, one can 
assume a coherent regime of multiple interactions only for not very large transverse momenta, which is restricted at 
mid rapidities by 

PT<^^. (15) 
mNRA 

At the RHIC energy, ^/s — 200 GeV, only quarks with up to several GeV transverse momentum can be produced 
coherently. 

The transverse momentum distribution of a parton after propagation through a nucleus is given by the square of 
the amplitude of multiple interactions. The quark impact parameters in these two amplitudes, direct and conjugated, 
are different. As a result, one can express the p^-distribution in terms of the eikonalized partial elastic qq dipole 
amphtude p7j . 

^^!Mll|^(6,p^,^) = -1^ / dMdV2e'^^-(---)l7?„(fi,f2)(l-Im/^^-(6,(fi -r^),/?)), (16) 
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where flj^{fi, is the density matrix which describes the impact parameter distribution of the quark in the incident 
hadron, 



nUrur2) = ^-^e-^^^'^+^'^^^^'-\ (17) 



where (fc|,) denotes the mean value of the parton primordial transverse momentum squared. The function Im/^ in 
the above denotes the partial amplitude of a qq dipole colliding with a nucleus at impact parameter b and can be 

fN 
Iqq 

impact parameter 6, 



written, in the eikonal form, in terms of the dipole elastic amplitude /i^ of a qq dipole colliding with a proton at 



Im/,^^(6,r,/3) = 1 - exp[- j d'slinf^^{s,r,(3)TA{b + s)], (18) 

where Tyi(6) = J dzpA{b,z) is the nuclear thickness function and pA{b,z) denotes the nuclear density at a impact 
parameter b and a longitudinal coordinate z. The fractional light cone momentum x of the target gluons is implicit 
in the above expression, and the parameter (3 in the qq dipole amplitude of the saturation model is taken to be 1/2. 

The integral over fl in Eq. (jl6p can be readily done. In order to compute the remaining integrals, we choose a 
coordinate in which the angle between the impact parameter b with s,r — r\—f2 and px are denoted by 9, 6 and 0, 
respectively. Then equation (|16p takes the following simple form. 



df3rdre^' 



(19) 



with the notation, 

1(6, r,S) ^ I d' s lmL%{s, f)TA {b + . (20) 



The result of the convolution of the nuclear thickness with the dipole amplitude will explicitly depend on the angle 
5 between the impact parameter b and the dipole transverse vector r. The azimuthal asymmetry resulting from a 
single quark passing through the nucleus is computed as a second order Fourier coefficient in a Fourier expansion of 
the azimuthal dependence of a single-particle spectra Eq. around the beam direction, 

r d^cos(20)^^^^Mziei 

^2bT, b) = da^(qA-,qX ) ■ (21) 

Substituting the expression in Eq. (fT9|) into Eq. ((2T|) , the integral over (j> can then be performed analytically by using 
the identity given in Eq. (jASp . getting 

/j" M /" rdrJ„lj,Tr)e-^'''-i<':'-'> 

where Jn{x) denotes the Bessel function. In the above expression the angle dependence 4> between the impact 
parameter 6 and the transverse momentum of the projectile quark pT disappeared and the azimuthal asymmetry is 
directly related to the dipole orientation with respect to impact parameter b through the angle 5. If one neglects the 
angular dependence of dipole cross section, then u| becomes identically zero regardless of a given nuclear profile and 
of dipole amplitude parametrization. 



V. V2 OF GLUONS 



One could calculate the cross section of high-pT gluon production in the same way as was done for quarks. Namely 
gluons can also experience multiple coherent interactions provided that the final pt is restricted by the condition 
Eq. (jlSp . One can neglect interaction with the spectator partons, like we did it for quarks in the previous section, 
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if the final px of the parton is much larger than its primordial transverse momentum. Otherwise, interaction with 
spectators is important since color screening is at work. Such an approximation is valid for valence quarks, which 
have a primordial momentum of the order of Aqcd, starting from px of several hundred MeV. 

Nevertheless, the dynamics of gluon radiation is more involved. There are many experimental evidences (l^ . [isj 
for a much higher primordial momentum of gluons compared to quarks. Therefore, one can neglect spectators only 
at pT of several GeV. In order to be able to work at smaller px one should include interaction with spectators, i.e. 
instead of "elastic" gluon scattering, GN — > GX, we need to consider bremsstrahlung subprocesses, GN 2GX, or 
qN qGX. The Born approximation for this processes includes three graphs (interactions with the initial and two 
final partons). The cross section of this reaction can be also expressed in terms of the dipole approach [3, [l^, and 
can be eikonalized on a nuclear target provided that the coherence length is sufficiently long. 

The condition of coherence is derived differently from the analysis that lead to Eq. ([U]). In this case Ic is simply 
the inverse longitudinal momentum transfer, 

^^=li' (23) 
where E is the initial parton energy, and M is the invariant mass of the two final partons, which reads, 

= -f—. (24) 
a[l — a) 

Here px is the relative transverse momentum of the final partons, a is the fractional light-cone momentum of one of 
the final partons, and parton masses are neglected. Since gluon radiation is dominated by small values of a ^ 1, 
Eqs. and lead to the same expression for the coherence length Eq. (fTi|) . where Eg should be replaced by the 
energy of the parton detected in the final state. 

In the long coherence length (LCL) regime, the transverse momentum spectra of gluon bremsstrahlung for a high 
energy gluon interacting with a nucleon N (nucleus A) target including the nonperturbative interactions of the radiated 
gluon reads [3, , 



da{GN{A) G1G2X) _ 1 



d^prd^b 2(27r) 

X Im [/3^J^) (6, fi, x) + /3^(^) (6, f2, x) - /3^(^) (&, (ri - rl), x)] , (25) 

where a — P+{Gi)/ P+{G) denotes the light-cone momentum fractional of the radiated gluon. The partial amplitude 
can be given in terms of the qq dipole amplitude, 

Imf^aib, r, x) = ^{Im/,^(&, f, x) + Im/,^(6, ar, x) + Im/,^(&, (1 ~ a)f, x)}, (26) 

where the factor 9/8 is the ratio of Casimir factors. Here the vectors r, ar and (l — a)r denote the two gluon transverse 
separations r(Gi) — r{G2), r{G) — r[G2) and r(G) — r(Gi), respectively. Notice that the parameter /3 that is present 
in the dipole saturation model of Eq. ([8|) corresponds here to the fraction of the total 3G momentum carried by the 
2G = G — G2 system, and is related to the light-cone momentum fractional of the radiated gluon a by 

/? = l-f- (27) 

We still have to specify the light-cone distribution hmction (^'gg) for GG Fock component fluctuations of the 
incoming gluon, which includes nonperturbative interactions of these gluons. This is given by. 



^'gg(^, a) 



2rl 



a{el ■ e)(e2 • r) + (1 - a)(e2 • e){el ■ r) - a{l - a){el ■ e2*)(e • r) , (28) 



where tq = 0.3 fm is the parameter characterizing the strength of the nonperturbative interaction, and which has been 
fitted to data on diffractive pp scattering [1^]. In Eq. P5|) the product of the wave functions is averaged over the 
initial gluon polarization, e, and summed over the final ones, ei_2- 

We consider here a case relevant for high energy gluon radiation with a — *■ 0. The azimuthal asymmetry of gluons 
V2 coming from a gluon- nucleon collision can be defined in a similar way to the given in Eq. (|2ip , although replacing 
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the particle spectra with the one in Eq. ((25)l . After some algebra one obtains, 

gN , ,N J-TT ^ d^pTil^b 



J,°"drj:^d6 cos(2<5)Im/3^(6,f)j 




^/(r,<5)| 


J,°^drj:^dS Im/iV(^-lj 


(1 - e-PT-oV2) j^(py^)e^ _ JobTr)e^/(r,<5)j 





(29) 

where the function /(r, (5) is defined as 

f{r,6) = dA de-—-^:^——^ ,,,, e ^. (30) 







(A2 +r2)2 -4(Arcos((5-6l))2 



In the case of a nuclear target Eq. ([25|) still holds, but the dipole amplitude on a nucleon target J^q should be 
replaced with the one on a nuclear target f^Q. The partial elastic amplitude f^Q for a colorless three- gluon system 
colliding with a nucleus can be written in terms of the partial amplitude f^Q of a three-gluon system colliding with a 
proton at impact parameter 6, 



Im/a^j (5, f, a;) = 2 1 1 - exp[- j d^s Ini/3^ (s, r , x)Ta (6 + s")] | , (31) 

where the 3G amplitude J^q is related to the qq dipole amplitude via Eq. ([26]). In a very similar fashion, one can 
obtain the gluons V2 in a gluon-nucleus collision (see Appendix A for a derivation), 

gA, ^ cos(20) '^-^'^,y;.f^^ ^ /-dr/:^d^cos(2J)v&^(pT,r,b,J) 

^ d<^ '^'J.pJ'f^'''' r dS ^niPT, r, 6, S) + 2ng{pr) ' ^''^ 

(1 - e-^-'-o/^) Ji(pTr)e4"^"^''"''^ + 2^Jo(pTr)e^"^"'''"''V(^, 5) 

(33) 



where we defined 




^N{PT,r,b,6) = 


(27r)2 , 


PT ^ 


*£,(PT,r,6,(5) = 


2(27r)2 


PT 


5(pt) = 


''f (1 






with the notation. 





Xg (&, = i y ^ W,,- (s, r)r^ (5 + , (34) 

The remaining integrals in Eqs. (|29l30l32l34p can be performed only numerically. From Eqs. (|29l32p one can observe 
that similar to Eq. (|22p . the angle dependence (j) between the impact parameter b and the transverse momentum of 
the projectile gluon px was replaced by the angle 5 between b and dipole vector r. As a consequence the azimuthal 
asymmetry of the radiated gluon is directly related to the orientation of the color dipole. 

VI. PIONS AZIMUTHAL ASYMMETRY IN pp AND pA COLLISIONS 

The invariant cross section for hadron production in pp collisions can be described, in the pQCD-improved parton 
model based on factorization [2l|, [2^ , by the expression 



= [ dx,dxJd^hTd^kJTf^/p{x^,Q^)gp{k,T,Q^)fJ/p{xJ,Q^)gp{kJT,Q^)K^{ij kl) '^Vfc(^fc'Q ) ^ 

ijkl 

(35) 
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where we sum over different species of participating partons, and fi/p{xi, Q^), fj/p(xj, Q^) are the parton distribution 
functions (PDF) of the coUiding protons, which depend on the hght-cone momentum fractions Xi, Xj and the hard 
scale Q. The function Dfi/k{zk, Q^) is the fragmentation function of parton k to the final hadron h with a momentum 
fraction z^. The cross section ^{ij — + kl) of the hard process, which is a function of Mandelstam variables, can be 
calculated perturbatively. 

The higher order perturbative corrections are taken into account via a _K-factor, and the primordial momentum 
distributions GpikT^Q"^), which are assumed to have the form 

yp[kT,Q) = ' 

with the mean values (fc|,), are taken to be independent of for the sake of simplicity. 



A. Short coherence length regime 

An important QCD prediction is the px power dependence of jet production, confirmed by data. Due to this property 
the multiple interactions of a large pt produced parton do not share equally the total transferred momentum pT (like 
it would be if the pt dependence of each collision were Gaussian). In fact, there is one collision with a large transverse 
momentum, close to pt, while the others are mainly soft interactions with small transferred momenta. Equation 
(I16p fully includes this dynamics, although it is really valid only for coherent multiple rescatterings, when multiple 
interaction amplitudes, rather than cross sections, are convoluted. What happens if one decreases the energy, or 
increases px, and eventually gets into a regime of short coherence length (SCL) (fT4|) ? According to the specific QCD 
dynamics described above, only the single high-pr collision becomes incoherent, while the other multiple soft collisions 
remain coherent. Thus, we arrive at a three step picture: soft multiple coherent interactions slightly increasing 
the transverse momentum of the parton, followed by a hard incoherent parton- nucleon collision with high px, and 
eventually multiple soft final state interactions of the produced parton leading to an additional broadening. A proper 
tool for calculations of soft multiple broadening is the dipole approach , which allows to use the phenomenology 
of soft interactions. At the same time, for the hard parton-nucleon collision we rely on the factorization based parton 
model, since at large pt Bjorken x of the target is too large for the dipole technique to be valid. 

Thus, the inclusive cross section of pA — s- hX can be written as 



dy(PpTCPh 



= TA{h) ^ J dXidXJd^hTd^kjTf^/pix^,Q'^) Gp{kiT,Q'^)fj/NiXj,Q'^)GpikjT, Q^) 

ijkl 

^K^{zj->kl)^I^l^^^^. (37) 

dt TTZk 



The parton distribution function of a bound nucleon, fj/N{xj , Q^) , is known to be modified by the nuclear environment, 
a phenomenon called EMC effect [IJ]. Notice that shadowing effects should not be considered, since we need the PDF 
of a single bound nucleon. Moreover, in the SCL regime the coherence length is too short for any shadowing effects 
to appear. At large Xj isotopic effect may be important for the target nucleon PDF, therefore we average over the 
nucleus, 

fj,N{x, Q2) = |/,/p(x, Q2) + - I) f,/n{x, Q"), (38) 

with atomic and charge numbers A and Z respectively. 

Initial/final state broadening of the projectile/ejectile partons is effectively taken into account via a modification 
of the primordial transverse momentum distribution, Gpikir) ^ Gpikix), where 

Gpik..) = (39) 

The /cT-distribution dN^"^^"^-^ [b) / d'^ kx j normalized to unity, is calculated using Eq. and the KST parametrization 
of the dipole cross section (with a Casimir factor 9/4 for gluons). Notice that in the SCL regime under consideration 
the process of broadening is dominated by soft coherent multiple interaction which have no relation to the hard scale 
imposed by the high-pj- process occurring incoherently. Therefore in this case the density matrix Eq. ([T7]) should 
have a rather small mean value of the primordial quark momentum (A:^) ~ ^qcz?- '^'^ simplify the calculations we 
assume that the initial and final partons are the same, so the total nuclear thickness TA{b) contributes to broadening. 
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B. Midrapidities at high energies 



At high energies and midrapidities the parton fractional momenta in the beam and target are small, xi ^ X2 ^ 
^VtI \fs <C 1, so hadron production is dominated by fragmentation of radiated gluons, and we can rely on the results of 
Section V. The cross section of hadron production in collisions at impact parameter 6 is then given by a convolution 
of the distribution function of the projectile gluon inside the proton with the gluon radiation cross section coming 
from GN collisions and also with the fragmentation function, 

/• , , , ^.^ dajGp ^ G,G,X) Dh/gA^,Q^) .... 
l^^^b = J ^^^^-/^(^-^ ) d^k,rd^b 7^ ■ (4°) 

To simplify the calculations we assume here that the projectile gluon has the same impact parameter relative to the 
target as the beam proton. At midrapidities we have Xg = 2kgT/\/s and we take = k^j,. The cross section of 
gluon radiation in the above expression can be obtained by the master Eq. This cross section reproduces well 

the measured pion cross section in pp collisions [20!|. In the LCL regime, a high pT parton propagating through 
the nucleus is freed by the multiple coherent interactions, and the Cronin effect may be conceived as color-filtering. 
In this regime the cross section of hadron production in pA collisions has the form, 

- ' dxgfG/p{xg,Q ) ,2, ,2, ^2 ' (41) 



dyd^PTd^b J " ' d^kgTd% 

where the cross section of gluon radiation in GA collisions can be obtained from Eqs. (|25I31[) . In Eqs. (|40I41|) . we will 
use the same PDFs and the fragmentation functions which are employed in Eq. (I35p . 



C. Azimuthal asymmetry of produced hadrons 

Notice that the asymptotic expressions (|40[ HT|) . supplemented with the gluon radiation cross-section given in 
Eq. (|25|) at a ^ 1, are only reliable at very long coherence lengths, which is certainly the case at LHC energies. At 
RHIC energies, for hadrons produced at midrapidities with moderate pt, we are in the transition region between the 
regimes of long and short coherence lengths. However, for peripheral collisions where V2 is not zero, we are in the 
LCL regime. 

The azimuthal asymmetry of hadrons produced in pp and pA collisions at impact parameter b is computed as a 
second order Fourier coefficients in a Fourier expansion of the azimuthal dependence of the transverse momentum 
spectra of the inclusive hadron production, Eqs. P7l40l4ip . around the beam direction. 



T 



vUPT,b)^ (.rj; ■ (42) 



dyd^P'j'd'^b 

In the SCL regime, the main element of the hadron production Eq. (|37p. which leads to an azimuthal asymmetry, is 
the angular dependence of the broadened projectile partons via Eq. (fT9|l . At the same time, in the LCL regime, the 
angle dependence of the gluon radiation cross section Eq. (|25p induces via Eqs. (|40l4ip an azimuthal asymmetry for 
the produced hadrons. In both cases, the introduction of the color-dipolc orientation is the key ingredient. 



VII. NUMERICAL RESULTS 



First, we show the V2 of quarks and gluons in parton-nucleus collisions introduced in sections IV and V. The only 
external input here is the nuclear profile. We take first the Woods-Saxon (WS) profile, with a nuclear radius Ra — 6.5 
fm and a surface thickness ^ = 0.54 fm. In our approach, the profile of nuclear density at the edge is a very important 
input, since the elliptic asymmetry stems from the rapid change of nuclear density at the edge. In order to show this 
more clearly, we also show the results for the hard sphere (HS) nuclear profile with a constant density distribution. 
Pa — Po^{Ra — f)i with the same nuclear radius as was taken for the WS profile. 

In Fig. ([1]), we show V2 of quarks at various impact parameters, within the saturation model I, and at the RHIC 
energy. A smaller primordial transverse momentum (k^) leads to a bigger broadening and more multi-scatterings, 
consequently the azimuthal asymmetry will be also bigger. 

The main source of azimuthal asymmetry in the amplitude p8l3ip is the interplay between multiple rescatering and 
the shape of the physical system. The key function which describes the effect of multiple interactions is the eikonal 
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FIG. 1: Right: The azimuthal anisotropy of quarks from quark-nucleus collisions, obtained from Eq. (I22|l . for various primordial 
transverse momentum squared {fcl^) at a fixed impact parameter b = 7. Left: The azimuthal anisotropy of quarks from quark- 
nucleus collisions at various impact parameters, with a fixed primordial transverse momentum squared {k^) — 1 GeV^. All 



curves are obtained for the RHIC energy ^/s 
profile. 



200 GeV, within the GBW dipole model with Woods-Saxon (WS) nuclear 



exponential and the color dipole amplitude. The information about the shape of the system is incorporated through 
a convolution of the impact parameter dependent partial elastic amplitude and the nuclear thickness function. The 
initial space-time asymmetry gets then translated into a momentum space anisotropy by the double Fourier transform 
in Eqs. p6l I25p . For more central collisions, the correlation between nuclear profile and dipole orientation is minimal. 
In fact if the nuclear thickness was constant, then the convolution between the nuclear profile and the dipole orientation 
would become trivial and there would be no azimuthal asymmetry. Therefore, the main source of azimuthal anisotropy 
is not present for central collisions where the nuclear density has only small variation. This can be seen in Figs. (|ll2p . 
where a pronounced elliptic anisotropy is observed for collisions with impact parameters close to the nuclear radius 
Ra, where the nuclear profile undergoes rapid changes. In Fig. ([2|), for comparison, we show the ^2 coming from the 
WS and the HS nuclear profile, within the GBW model. The u| of quark with the HS nuclear profile can be twice 
bigger than the one with WS nuclear profile. This indicates that the nuclear density profile is an important input. 

In Fig. ([2]), we show the of quarks in pA collisions within the GBW model and the improved Born approximation, 
the 2g model. For the 2g model, we show for two different effective gluon masses nig = 140 and 650 MeV. There 
are growing evidences suggesting that the effective gluons mass should be bigger than the confinement scale. For 
instance, a small gluon correlation radius 0.35 fm is a result of lattice calculations [1^, and it is also predicted by the 
instanton model [26| . Experimental data suggests an enhanced intrinsic motion of gluons in light hadrons compared 
to the inverse hadronic radius [l9[. And also the smallness of the triple-Pomeron coupling can be only explained by 
a enhanced gluon interactions due to nonperturbative effects p^ . 

It is obvious from Fig. Q that the azimuthal asymmetry of quarks is very sensitive to the effective gluon mass 
and the higher order corrections. The azimuthal asymmetry within the saturation model which includes higher order 
radiation corrections is significantly smaller than the one for the 2g model. Therefore, the inclusions of the higher-order 
radiation corrections or taking a bigger effective gluon mass reduces the azimuthal asymmetry. A quark propagating 
through a nucleus interacts by gluon exchange with different nucleons located at different azimuthal angles relative 
to the quark trajectory, and their contributions to V2 tends to cancel each other, reducing the azimuthal symmetry. 

In Figs. pi4p . we show the azimuthal asymmetry of gluon radiation for gluon- nucleon and gluon-nucleus collisions, 
given in Eqs. (|29l32p at various impact parameters. The azimuthal asymmetry of gluons in both GN and GA has a 
rather similar trend; it can be negative at small pr for peripheral collisions and becomes positive at higher pT- Again, 
for more peripheral collisions where the color dipole orientation becomes more important, the azimuthal asymmetry 
is bigger. In Fig. on the left panel, we also show the gluonic V2 in GA collisions at a fixed impact parameter b = 7 
fm, within various color dipole models. Notice that for the case of gluon bremsstrahlung we have already included 
some nonperturbative effects of gluon interactions through the GG distribution function Eq. (j28p . via the parameter 
rg = 0.3 fm, which simulates the strength of the non-perturbative gluon interactions contained in the diffractive 
pp scattering data. Therefore, the azimuthal asymmetries of gluons within the saturation dipole model and the 
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FIG. 2: The azimuthal anisotropy of quarks coming from quark-nucleus scatterings, for both the Woods-Saxon (WS) and 
hard sphere (HS) profiles, within the GBW dipole model and the improved Born approximation model, at the RHIC energy 
"v/s = 200 GeV. For the 2g model, we show the u| for two different effective gluon masses mg = 140 and 650 MeV . 




FIG. 3: The azimuthal anisotropy of gluons from gluon-nucleon collisions, given by Eq. (|29|) . at various impact parameters, 
within the GBW dipole model, and with the WS nuclear profile. 



improved Born approximation with effective gluon mass mg — 650 MeV are rather similar. 

Next, we calculate the azimuthal asymmetry of produced pions in and collisions, within two different 
schemes, the standard pQCD factorization Eqs. (I35I37|) (short-coherence) and light-cone factorization Eqs. ()40I41|) 
(long-coherence) . 

We first concentrate on the short-coherence scheme. For the parton distribution fi/p{x, Q^) used in our calculations, 
we employ the MRST (2006 NNLO) set [l^. For the fragmentation function D{z,Q^), we use the parametrization 
given by Kretzer [2^, with NLO corrections. In our calculation we will take (fc^) = 3 GeV^. This value corresponds to 
an average transverse momentum of the parton fcy ~ 1-5 GeV, which coincides with the analysis of fc^ for photon and 
pion production given in Refs. [1^, [sO]. Notice also that different value of has been used in different approaches 
[3l|. For the scale Q of the hard process in Eqs. p5l37p . we choose Q = pt/{1-2zc). We take the if-factor K 2, 
which gives a good approximation of the NLO order contribution in the px region of interest. With this setup, we 
are able to describe pion production data for pp collisions from the SPS to RHIC energies within 35% discrepancy, 
see Fig. ([5]) right panel. 
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FIG. 4: Right: The azimuthal anisotropy of gluons from gluon- nucleus collisions, given in Eq. (|32|l . at various impact parameters, 
within the GBW dipole model. Left: w| at a fixed impact parameter 6 = 7 fm, within various dipole models. We used for all 
curves the WS nuclear profile. 
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FIG. 5: Right: the cross section of p -j- p ^ 7r° -|- X at RHIC for y/s = 200 GeV, and that of p -j- p ^ 7r+ -j- X at the SPS fixed 
target experiment for Eiab = 300, 400 and 800 GeV. Data for RHIC and SPS are from Ref. [33] and Refs. [33l. [34|. respectively. 
Left: The azimuthal anisotropy of the produced pions in pA collisions, in the short-coherence regime, at energy Eiah = 400 
GeV, for the Woods-Saxon (WS) nuclear profile, and within the KST dipole model. 



In this way, all free phenomenological parameters for hadron production in yA collisions, Eq. (|37p . are already 
fixed in reactions different from yA collisions. We recall that the azimuthal asymmetry originates here from the 
rapid change of the nuclear density, and therefore only the peripheral tail of the nuclear profile contributes to the 
azimuthal asymmetry. At the very periphery, the nuclear parton distribution is unchanged compared to a free nucleon. 
Nevertheless, we have numerically verified that changes less than 20% with various PDF parametrizations for impact 
parameters bigger than 5 > 6 fm. In Fig. ((5|) we show the azimuthal asymmetry of the produced pions in pAu collision 
at midrapidity, for the fixed target energy Eiab = 400 GeV obtained in the SCL scheme, Eq. ([37]). 

In Fig. ([S]), right panel, we show the azimuthal asymmetry of produced pions in pp collisions at the RHIC energy 
= 200 GeV. A pronounced positive azimuthal asymmetry is observed, at an impact parameter around half of the 
proton size h ~ 0.5 fm. The source of the azimuthal asymmetry of pions in collisions is the angle dependence of 
the gluon bremsstrahlung cross section in Eqs. via the color-dipole orientation. The azimuthal asymmetry of the 
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FIG. 6: Right: The azimuthal anisotropy of produced pions V2 in pp collision at the RHIC energy — 200 GeV, at various 
impact parameters. Left: The azimuthal anisotropy of produced pions in pAu collision with the WS nuclear profile at the 
RHIC energy. We used for all curves the GBW dipole model. 



fragmented pions in pp collisions is sizable at a pt where the color-dipole size becomes compatible with the impact 
parameter and consequently the color-dipole orientation becomes important, sec Eq. (|5]). 

At the LCL hmit, pions are entirely produced by gluons. It may be then puzzling that the pion's V2 is positive 
while gluons have negative V2 in the same range of transverse momentum. Nevertheless, one should notice that the 
gluons transverse momentum kgT is related to the transverse momentum of the fragmented pions via kgx = Pt/z. 
Therefore, high kgx gluons with positive u| are responsible for the produced pions at moderate px, see Fig. ([3]). 

In Fig. ^ left panel, we show the of the produced pions in pAu collisions, obtained in the LCL scheme. The 
azimuthal asymmetry increases from zero at an impact parameter around nuclear radius to a maximum value at 
an impact parameter around 6 ~ 7 fm, and then decreases again for a more peripheral collision. 



VIII. SUMMARY AND CONCLUSIONS 



This paper is the first attempt to calculate the azimuthal asymmetry V2 of produced hadrons in pp and pA collisions. 
We generalized the hadron production scheme within the pQCD parton model and the light-cone QCD-dipole formal- 
ism by the inclusion of the color-dipole orientation. This is the key ingredient which leads to the azimuthal asymmetry 
both in pp and pA collisions. We introduced the color-dipole orientation into the improved Born approximation and 
within the saturation model of Golec-Biernat and Wiisthoff, which satisfies available DIS data. 

In the short-coherence regime, we used the improved pQCD parton model. The salient part of hadron production 
in pA collisions at the short-coherence regime is the broadening of the projectile partons going through nucleus. This 
broadening is calculated here in a free-parameter formalism which uses the color-dipole approach, where the angle 
dependence of the broadened projectile partons via the color-dipole orientation leads to the azimuthal asymmetry of 
the produced hadron. In the long-coherence regime, we used the light-cone dipole formalism in the rest frame of the 
target, which is valid at small X2- The main source of the azimuthal asymmetry in the long-coherence regime originates 
from the angle dependence of the radiated gluons cross section. Our results show that the azimuthal asymmetry of 
pions in both pp and pA collisions is rather small. This indicates that the contribution of the initial state effects, which 
is present in cold nuclear matter in the observed azimuthal asymmetry V2 of the produced hadrons in AA collisions 
at RHIC, is very small. 

We have also systematically studied the azimuthal asymmetry of partons in parton-nucleus collisions. The azimuthal 
asymmetry of partons in partons-nucleus collisions is very sensitive to the effective gluon mass and the higher-order 
radiation corrections. We found that the azimuthal asymmetry of gluons at small px can be negative in Gp and 
GA collisions, and at higher px changes sign and becomes positive. This is in contrast with the ?j| of quarks in qA 
collisions, where it is positive. 

The technique presented in this paper can be also used to study the azimuthal asymmetry in DIS and in the 
production of dileptons. We plan to report on some of these problems in the near future. 
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FIG. 7: The angles defined in Eqs. (|A4IAlip . 
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APPENDIX A 



In this section we illustrate a derivation of Eq. At a = 0, up to a constant normalization N, the averaged 

product of the light-cone wave function Eq. ([28]) reads, 

^^^-^^^ (Al) 



2rg 



2 2 



^'QQ(ri,a)*GG(r2,a) = TVexp 
The azimuthal asymmetry is defined as, 

where the gluon radiation cross section is defined in Eq. ([25]). The Fourier-integral is inconvenient for numerical 
calculation, but we can perform some of integral analytically in order to make the numerical task manageable. Let's 
define 

Ig {b, r,S) ^ ^1 <fs Im/^^ (s, r)TA {b + s), (A3) 

where the factor | comes from the definition Eq. (j26p at a = 0, and S denotes the angle between the impact parameter 

b and the dipole vector r, see Fig. [71 

We first perform the calculation for the first term in the numerator and the denominator /^vi, Idi in Eq. (jA2p . 



'1'2 



= -2 / d(j)dridr2dddd2 I ^ M cos(5 - 6*2)6 ^"0 e ^"0. (A4) 

The angles in the second Eq. (|A4|) are defined in Fig. [T] We can immediately perform the integral over 62 in the above 
equation by using the following identity 



-t-00 

^iAcos{(p) 



e 

n— — 00 



J2 ^"J„(A)e™^ (A5) 
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where J„(a;) denotes the Bessel function. We have 



+ OC 



i— — oo 

— — 27ri Ji(pt'"2) cos(0 — 5). (A6) 
Having plugged the above equation into Eq. (jA4[) . we can then perform the integral over (j), 



(Jiiprri) - J3(PTfi))cos(2(5) 
2Ji{pTri) 



(A7) 

Now we can simplify Eq. (jA4p by using ICi and IC2 given in Eqs. (jA6IA7p . 

'^^1^ - -47r2 /d^dril^^'^i^^^'^i)^-;^^^^^^^^^^ ldr2MpTr2)e^ 



(27r)2 



The remaining integrals in the above expression can be only done numerically. Now we calculate the second term in 

'4 



the denominator and numerator of Eq. (|A2[) 



(A9) 



We first change the variables in the above integrals by: 



n = r2 = (AlO) 



Therefore, we obtain, 



ID2J ij V 1 / (A2 +r2)2 -4(A.r)2^ 



cos 



drdAd6dH02 [ -f^^) Jf .^.2 - (AH) 

V 1 / (A2+r2)2-4(Arcos((5-6l2))2 ^ ^ 

The angles in the above equation are defined in Fig. [71 Notice that the factor i in the above expression is the Jacobi 
determinant. The integral over (f) can be immediately performed using the identity Eq. (|A5p . 

-^Af2 7 -^£)2 defined in Eq. (jAll[) can be then written in the following form 

'^j^J^ = 2nJ drdS {^'^'^^j^ll'^f^^^) e^e-^-('^^''^fir,S). (A13) 
where the function /(r, 5) is defined, 

f{r,S)^ dAd02j-r^ ^-^ -T-c-^ 77 TTTT^e^. (A14) 

^ ' ^ J (A2 + r2)2 - 4(Arcos((5 - 6'2))2 ^ ^ 
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The remaining multi-integrals in the above equation can be only carried out numerically. The last term in the 
denominator of i;| is independent of impact parameter 5, 



lD3 



^ ri-r2 



' 1'2 



= 271 dndr2d0id02cos{ei)cos{e2)e ' " ' ' (A15) 

Using the following integral, 

J dOe'P^''' cos(0i) = 2TriJi{pTr). (A16) 



we obtain 

2 



Id3 = (2^)' drMpTr)e^^ = ^ (l - e-P^o/^Y . (A17) 



To put all equations together, we can write defined in Eq. (jA2p as f f = ^2 in a factorized form, 
„P = (i^e-p'^^'-'/^y -^^^ (l~e-P>'/^) f drJi{pTr)e'^ fo{b,r) + {2n) f drjQ{pTr)e^ fi{r,b), 



Pt ^ ^ PT 



Pt 



(1 - e-P-'-o/^) / dr{J^{pTr)-UpTr))e^ f2{b,T)^{2iT) / drJ2{pTr)e^ hM, 



(A18) 



where function /(r, 5) is defined in Eq. (jA14[) and /o_3(6, r) in the above equation are defined as follows 



/o(&,r) - 


j d5e-^'='^''-'''^\ 


hih.r) = 


j d6e-^<=^^''-''-^'^fir, S), 


f2{b,r) = 


j d5cos{25)e~^'^^^''-'^\ 




j d6cosi26)e-^'^'-''^'-'^^f{ 



(A19) 



In a very similar fashion one can also obtain Eq. ([29] 
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